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1  Summary 


This  project  has  investigated  the  application  of  wavelet  methods  in  spread  spectrum  communi¬ 
cations.  Use  of  the  wavelet  transform  as  an  alternative  to  the  short-time  Fourier  transform  and 
quadratic  time-frequency  transforms  (e.g.,  the  Wigner-Ville  distribution)  for  measuring  critical  pa¬ 
rameters  of  intercepted  slow  frequency  hopping  communication  signals  has  been  explored.  The 
wavelet  transform  was  found  to  not  offer  significant  advantages  over  Fourier  methods  in  this  ap¬ 
plication  because  the  sinusoidal  analysis  provided  by  the  Fourier  transform  is  well  matched  to  the 
piecewise-sinusoidal  carriers  used  in  frequency  hopping  signals.  However,  use  of  the  wavelet  trans¬ 
form  with  sinusoidal  wavelets  offered  an  efficient  mulliresolntion  framework  for  precise  measurement 
of  hopping  limes  and  frequencies.  A  hybrid  wavelet- Fourier  algorithm  was  developed  to  exploit  this 
feature. 

A  new  approach  to  multiple-access  spread  spectrum  communications  was  also  developed  under 
this  project.  This  approach  employs  crthogoral  wavelc'r  as  symbols  for  encoding  digital  information 
signals.  Several  users  can  transmit  simultaneously  in  the  same  medium  using  channels  defined  by 
wavelet  scale  in  much  the  same  way  that  channels  are  defined  by  frequency  in  traditional  frequency 
division  multiple  access.  With  appropriate  choices  of  wavelet  symbols,  this  approach  has  been  shown 
to  provide  spectral  spreading  that  makes  the  transmitted  signals  robust  with  respect  to  detection, 
interception,  and  jamming  by  typical  methods. 

2  Research  Objectives 

The  main  objectives  for  the  research  under  this  effort  are  summarized  as  follows: 

1.  Investigate  the  utility  of  the  wavelet  transform  in  estimating  the  spreading  function  for  inter¬ 
cepted  slow  frequency  hopping  spread  spectrum  communication  signals.  Identify  advantages 
and  drawbacks  of  the  wavelet  transform  compared  to  the  short-time  Fourier  transform  in  this 
application. 

2.  Develop  wavelet-based  approaches  to  multi-user  spread  spectrum  communications.  Consider 
the  design  of  wavelets  that  provide  robustness  of  the  transmitted  signal  to  detection,  intercep¬ 
tion,  and  jamming. 

The  results  obtained  in  pursuit  of  these  objectives  are  described  in  the  following  two  sections  of 
this  report'.  For  convenience  in  exposition,  objective  2  is  discussed  first. 
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3  Research  Results:  Wavelet-Based  Spread  Spectrum 


The  purpose  of  this  portion  of  the  research  effort,  was  to  introduce  a  new  method  of  coding  digital 
communication  signals  on  orthogonal  symbols  in  such  a  way  that  several  users  can  communicate 
simultaneously  using  the  same  transmission  medium.  Such  methods  are  known  in  the  communica¬ 
tions  literature  as  code  division  multiple  access  (CDMA)  schemes1.  The  two  most  popular  CDMA 
schemes  arc  frequency  division  multiple  access  (FDMA)  and  time  division  multiple  access  (TDMA). 
In  FDMA,  orthogonality  of  the  various  users’  messages  is  achieved  by  separating  them  into  essen¬ 
tially  disjoint  frequency  bands  while  orthogonality  of  TDMA  messages  is  obtained  by  constraining 
them  to  occupy  disjoint  time  intervals. 

Though  FDMA  and  TDMA  offer  advantages  in  simplicity,  they  are  obviously  not  the  only  ways 
to  achieve  orthogonality  of  transmitted  message  information.  Other  CDMA  schemes  were  introduced 
originally  in  the  context  of  military  communications,  where  they  were  used  to  provide  robustness 
with  respect  to  detection,  interception  by  unintended  receivers,  and  interference  -  including  deliber¬ 
ate  interference  (  jamming)  [20].  More  recently  the  value  of  certain  CDMA  approaches  in  commercial 
applications,  including  cellular  telephony,  has  been  pointed  out  and  civilian  CDMA  has  become  a 
topic  of  interest  within  the  communications  research  community. 

This  research  has  developed  a  CDMA  scheme  called  scale-division  multiple  access  (SDMA)  that  is 
based  on  the  theory  of  wavelets.  SDMA  offers  much  of  the  simplicity  of  FDMA  while  retaining  many 
of  the  desirable  features  of  more  complicated  CDMA  schemes,  including  robustness  to  detection, 
interception,  and  jamming.  In  FDMA,  bits  are  encoded  on  orthogonal  symbols  obtained  by  lime 
shifting  and  frequency  shifting  a  single  “mother”  symbol  -  a  truncated  sine  wave.  SDMA  encodes 
bits  on  orthogonal  symbols  obtained  by  lime  shifting  and  dilating  a  single  mother  wavelet  symbol, 
which  may  be  chosen  front  a  fairly  broad  class  of  wavelets.  It  is  the  separation  of  the  various  users’ 
messages  into  orthogonal  channels  at  different,  dilations  or  "scales”  that  suggests  the  name  “scale 
division  multiple  access”. 

Closely  associated  with  the  theory  of  wavelets  is  the  concept  of  a  frame,  which  generalizes  the 
notion  of  an  orthonormal  basis.  The  principle  of  SDMA  coding  has  been  shown  to  extend  from 
orthonormal  wavelet  bases  to  wavelet  frames. 

'CDMA  is  interpreted  in  *  broad  sense  here,  i  t.  as  coding  on  orthogonal  oi  nearly  orthogonal  symbols 
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3.1  Mathematical  preliminaries 

This  short  section  reviews  some  mathematical  concepts  and  introduces  some  notation  and  conven¬ 
tions  that  will  be  used  in  explaining  SDMA  and  in  the  subsequent  discussion  on  demodulation  of 
intercepted  frequency  hopping  signals. 

3.1.1  Lebesgue  spaces 

Unless  otherwise  noted,  the  expression 

Jsf(i)dt 

will  denote  the  Lebesgue  integral  of  the  Lebesgue  measurable  function  /  :  K  — >  C  over  the  Lebesgue 
measurable  set  S  c  8,  For  1  <  p  <  oo,  the  p  norm  of  the  Lebesgue  measurable  function  /  :  R  — *  C 
is  defined  as 

ii/iip  =  [JK\w)  i^y7*- 

If  the  integrals  of  two  Lebesgue  measurable  functions  /  and  g  on  every  measurable  set  are  equal, 
then  {  and  g  are  said  to  be  Lebesgue  equivalent.  T  his  condition  is  identical  to  the  requirement  that 
the  Lebesgue  measure  of  the  set  {t  \  f{t)  ^  g(t)}  is  zero. 

1  or  p  >  1 ,  the  ],ebesgue  space  Lv  ^  Lp(7.)  consists  of  equivalence  classes  of  Lebesgue  equivalent 
functions  /  for  which  ||/||p  <  oo.  With  addition  and  scalar  multiplication  defined  in  the  obvious 
ways.  Lp  is  a  complex  Banach  space  for  all  p.  Furthermore,  I.7  is  a  complex  Hilbert  space  with 
inner  product 

(/.<?)  =  J  f{i)g'(t)dt 

One  often  writes  “/  6  ip"  to  mean  that  the  equivalence  class  of  /  is  in  Lp. 

The  discrete  Lebesgue  space  lr  =  CP(Z)  consists  of  complex  sequences  c  =  { c/t  |  k  e  Z)  for  which 

I|C||„=  /P  <  OC 

\k(7.  J 

As  with  Lp ,  Cp  is  a  complex  Banach  space  for  ali  p  and  (7  is  a  complex  Hilbert  space  with  inner 
product 

(c,d)  =  YjCkdk 
ktl 

Lebesgue  spaces  and  norms  over  S’*  and  Zn  are  defined  analogously. 

Since  this  report  is  primarily  concerned  with  analysis  in  J.2,  ||/||  will  be  used  to  denote  ||/||2 
when  there  is  no  danger  of  confusion. 
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3.1.2  Fourier  transform 


The  Fourier  transform  of  /  e  ( L 1  n  L2)  is  the  funtlion  T{/}  =  /  C  with  values 

f(»)  =  I  Me-'**  dt  (1) 

JTL 

Lebesgue  equivalent,  functions  clearly  have  the  same  Fourier  transform.  Hence  the  operator  J  is  well 
defined  on  O  ft  L2.  It  is  clearly  iinear. 

Since  L1  0  L2  is  dense  in  L 2,  the  definition  of  7  extends  naturally  to  a  mapping  of  L 2  onto  L2 . 
The  particulars  of  this  extension  are  given  by  Planchcrcl’s  theorem  [23].  One  consequence  of  this 
extension  is  that,  for  functions  /  and  q  in  L2, 

<A$)  =  2»</.s) 

which  obviously  implies  ||/||2  =  27t||/||2. 

If  /  €  L2,  the  inverse  Fourier  transform  of  /  is  the  function  g  :  R  —  C  defined  by 


The  function  g  is  Lebesgue  equivalent  to  /.  Together  with  the  preceding  notes,  this  implies  that 
the  Fourier  transform  defined  by  (1)  extends  to  an  isomorphism  of  L2  onto  L2  which  is  isometric 
except  for  the  constant  factor  \/2x. 


3.1.3  Bandlimited  and  timelimited 

In  electrical  engineering,  it  is  common  to  regard  L2  functions  as  “signals”  which  are  functions  of 
time.  Then  Fourier  transform  variable  u>  is  thus  known  as  (angular)  frequency. 

A  function  /  €  L2  is  bandlimited  if  there  is  an  fl  e  1  so  that  /(u>)  =  0  for  u>  i  (-Sl.fi).  On  the 
other  hand,  a  function  /  el2  is  timelimited  if  there  is  a  T  6  E  and  a  function  g  that  is  Lebesgue 
equivalent  to  /  such  that  g(t)  =  0  for  i  t  (-7,7  ). 


3.1.4  Orthonormal  bases  of  L2 


If  J  is  an  index  set,  a  set  H  =  |  j  €  J)  of  elements  in  L2  is  orthonormal  if 


for  all  j  and  k  in  J . 


{<!>},  <t>k) 


1  j  -  k 

0  j*k 
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Figure  3:  An  example  of  encoding  information  on  the  Haar  wavelet  at  different  scales. 

robustness  with  respect  to  detection,  interference,  and  narrowband  jamming.  Design  of  wavelets  for 
SDMA  applications  is  discussed  in  later  sections. 

3.3.2  Channel  capacity  and  synchronization 

It  is  evident  from  Figure  3  that  transmitters  using  finer  scales  (lower  m)  will  be  able  to  transmit 
more  bits  per  unit  time  than  those  using  coarser  scales.  J'his  may  be  a  feature  in  some  applications; 
where  it  is  not,  time-division  within  the  high-capacity  channels  could  be  used  to  equalize  capacities. 

Expressions  (d)  and  (5)  assume  that  the  transmitters  are  perfectly  synchronized  and  the  prop¬ 
agation  delays  from  each  transmitter  to  the  receiver  are  identical.  Orthogonality  of  the  channels 
is  compromised  if  the  received  signals  are  not  perfectly  synchronized,  as  is  evident  in  the  example 
of  Figure  3.  This  problem  is  not  restricted  to  SDMA.  Rather,  it  is  inherent  in  any  system  using 
time-limited  symbols  [5].  Design  of  bandlimited  wavelets  for  SDMA  that  achieve  small  channel 
crosstalk  under  asynchronous  operation  is  discussed  later  in  this  report. 
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3.3.3  Coding  on  frames 


The  idea  of  coding  on  orthogonal  sets  can  he  extended  io  coding  on  frames,  t  he  extension  is  most  well 
suited  for  transmission  of  large  data  blocks.  Let  4>  =  {d>7  |  j  £  ./}  be  a  frame  and  let  =  jd>,  |  j  € 
be  its  dual.  If  b  =  |  j  6  J)  is  a  collection  of  scalars  corresponding  to  one  information  sequence 

as  in  part  A  of  this  section,  the  transmitter  forms  an  analog  signal  s  via  the  adjoint  frame  operator 
T'  operating  on  b 

-t*b  =  E 

,ieJ 

Without  noise,  the  receiver  obtains  a  corresponding  sequence  3  =  {3j  |  j  €  J)  by  applying  the  frame 
operator  T  to  the  received  signal  s. 


3  =  Ts  or  3j  =  (s,d>7) 

The  map  7T*  :  b  -*  3  is  not  generally  injective,  so  that,  an  equivalence  class  Bg  of  scalar  sequences 
is  mapped  to  the  same  3 ■  3  is  called  the  typical  element  of  Bg. 

If  different  information  sequence  corresponds  to  different  equivalent  class,  then  the  information 
sequence  is  one  to  one  corresponding  to  the  typical  element.  In  another  word,  the  receiver  can 
recover  the  information. 

In  1  he  presence  of  additive  noise  z ,  the  signal  at  receiver  becomes  s  +  z  and  the  demodulation 
operation  yields 

/'(#)  =  I{$  +  z)  =  3  +  Tz 

which  generally  doeb  not  correspond  to  the  equivalence  class  of  s.  To  combat  this,  a  superclass  can 
be  constructed  which  is  the  union  of  the  equivalence  classes  whose  typical  elements  are  near  in  the 
sense  of  the  norm  (note  T  is  continuous).  This  partitioning  represents  the  decision  mechanism. 

3.4  Spectrum  analysis  of  SDMA 

The  purpose  of  this  section  is  to  describe  explicitly  how  the  smoothed  power  spectral  density  of 
an  SDMA  signal  depends  on  the  mother  wavelet  and  the  statistical  behavior  of  the  information 
sequence.  A  procedure  for  pre-shaping  the  autocorrelation  function  of  the  information  sequence  is 
also  described. 
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3.4.1  Power  spectrum 

Beginning  with  the  form 

v{l)  =  Kc[s(«)c-''] 

which  relates  the  bandpass  signal  i>  to  the  wide-sonse  stationary  lowpass  signal  s.  the  autocorrelation 
'unction  of  v  may  be  expressed  as 


Ow(t)  =  Re[^.,(r)e,“'t7i 


where  oS3  is  the  autocorrelation  function  of  the  equivalent  lowpass  signal  s.  The  Fourier  transform 
of  <Pvt.  yields  the  desired  expression  for  the  power  spectral  density  4>vl  in  the  form 

-  u’c)  +  ~  w'c)] 

wlieic  is  the  power  density  spectrum  of  s.  Thus  it  sufTices  to  determine  the  autororrclation 
function  and  the  power  spectral  density  of  the  equivalent  lowpass  signal  s. 

First  consider  the  digital  modulation  methods  for  which  3  is  represented  in  the  general  form 


*(')  =  £<«(<) 
m,n 

where  b "  represents  the  n1*1  real-  or  complex-valued  information  symbol  in  t.e  with  channel  (see 
Figure  4)  and  IV  is  a  mother  wavelet.  The  set  {W"  |  2”m/?VF(^  -  n);n,m  e  1}  constitutes  an 
orthonorrnal  wavelet  set. 

The  autocorrelation  function  of  s  is 


4>(l  +  T\i)  = 


±E[a*(fM'4-r)] 

m'.n'.tn.n 


t  +  T 

2  m' 


(6) 
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Assume  that  the  sequence  of  information  symbols  is  wide  sense  stationary  with  autocorrelation 
function 

Olt(nr,v)=1-£[bX,b'^m}  (7) 

Then  (6)  can  be  expressed  as 

tf(f  +  r:<)=  E  0..X-  m;n'-n)2'l”*,+m»/2H"(^r  -n)W(~  -  n')  (8) 

Imposing  the  requirement  <p„(m;  n)  —  0  for  m  ^  0  (i.e.,  that  the  information  streams  from  different 
users  arc  uncorrelated)  simplifies  the  analysis.  Then  (8)  becomes: 


+  r;  1) 


=  E  2-^(0:  ’( ^r-~  ^  ) 

=  E  2-m6(0~.  n')Y  U"(  — 


The  summation  in  (Q),  namely, 


(9) 


Ew'*( 


f  -  2’"n  f  -  r  -  ?'”n'  -  2  nn  ^ 
2m  2’" 


is  periodic  in  the  t  variable  with  period  2"‘  in  the  channel  m.  But  0[t  +  r;t)  is  generally  not  periodic 
{if  it  is,  the  process  is  called  cycloslotwnary).  Such  a  process  having  multiple  periodicities  is  called 
a  pseudo-cyclostnhonary  process. 

In  order  to  compute  the  power  spectral  density  of  a  pseudo  cyclostationary  process,  the  depen¬ 
dence  of  d>,, (t  -f  r;t)  on  the  t  variable  must  be  eliminated.  This  can  be  accomplished  simply  by 
averaging  d>,,(t  -f  r\t)  over  a  period  within  the  different  channels.  1'hus 


=  E  <M0;n')2-m2- 

n',m,n 

=  E  0„fO;rr,)2  m2-m 

n',m 

=  E  <t>t,(0\n')2-m2-m 

n',m 

=  Y  </>,.(0;rT)2-m2-m 

n\m 


t*  i*  i  -  2"‘w 
J-7mJ 2  2 

?/- 
?/- 


i  -  2"‘n . t  +  r  -  2'7V  -  2^n 


2” 


)<*' 


2"7J  w-(  IzYJi )W/(  ^ 

jm/j  2’”  2m 

rn H-r'n  £  t  +  r  -  2mr?' 

7tn  l7-7nin  2m 

nw  ' 


(10) 
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The  Fourier  transform  of  the  relation  in  (10)  yields  the  power  spectral  density  of  *  in  the  form 


=  -F(<Mr)) 

=  £  <M0;  n')2~2rne~ '*2"‘ n'22rr‘\ l'V(2’r‘Lj ) | 2 

n\m 

-  ^  d„(0;n')e-“'2'"r‘'|Vi-(2mw)|2  (11) 

n#  ,»n 

where  IF  is  the  Fourier  transform  of  IF.  Defining 

*.,M  =  £fcf(0; 

u' 

t he  result  in  (11)  can  be  expressed  as 

4>,,M  =  ]>>„(2m^)|lV'(2’V)l2  (12) 

m 

Defining  5(u)  -  $„(u;)|M"’(u;}j2  allows  Kquatiori  (12)  to  be  further  simplified  to 

4>„U>)  =  £S(2"'w)  (13> 

m 

This  result  illustrates  the  dependence  of  the  power  spectral  density  of  s  on  the  spectral  characteristics 
of  the  mother  wavelet  IF  and  the  information  sequence  {6^}.  That  is.  the  spectral  characteristics 
of  ,,  ca(,  ho  controlled  by  design  of  the  mother  wavelet  IF  and  by  design  of  the  autocorrelation 
characteristics  of  the  information  sequence. 

3.4.2  Autocorrelation 

Since  the  spectrum  of  SDMA  signal  depends  on  the  autocorrelation  of  information  sequence,  it  is 
interesting  to  consider  shaping  of  the  autocorrelation. 

Cbim  3.1  Suppose  a  random  sequence 

. . .  h_ 2  6_i  6q  !>i  ... 

is  rtdr-seri.se  stationary  with 


{<p(0)  when  m  =  n 

i 

0  otherwise 


I<"  i\nl  rat  ms  A'o,  A-|, k2 . fc,v_ define  a  new  sequence 

...  r/_  2  7-1  <7o  7l  72  ... 

by 

q,  =  kob,  +  frift,+i  +  fc26t+2  +  . . .  +  A:ri_16t-.;./vr_1 

7V,en  [ft,-}  cart  6e  recovered  from  {7,}  6y  knowing  consecutive  initial  N -lb, ’s  and  the  autocorrelation 
{f>q}  of  the  sequence  {7,}  is  given  by 


*,[0] 

fro  kj  fr-2  ...  i*/v_ 2  k^-\ 

fro 

*,[1] 

k\  ^2  kj  ...  kfj- 1  0 

fr'l 

k 2  fr'3  k4  0  0 

fr"2 

<M*  -  1] 

fr  A/  _i  0  0  .  ... 

fr/V-1 

(M) 


and 


<t>,[k]  =  <fy[-fr] 

•  0,,(fr)  =  0  for  k  >  N 

The  proof  is  omitted  here. 

This  laim  gives  a  method  to  design  sequences  with  desired  autocorrelations  from  ;,}  by  select¬ 
ing  differ ont  k0,. . .,  fr^-i-  The  matrix  equation  (14)  can  also  be  written  as 


fro 

frl 

fr'2 

-II 

II 

fr’N- 1 

l 

fro 

frl 

fr'2  • 

•  fr/V-2 

fr'N-l 

-1 

<M0) 

frl 

frj 

kj  ■ 

. .  frN-| 

0 

**[1] 

frj 

fra 

k4 

..  0 

0 

^[2] 

frN-l 

0 

0  . 

*#[at  - 1) . 

(15) 


So,  conversely,  if  <£,[0],  <£?[1], ...  <f>,[N  -  1]  are  given,  fr0,  fri ,  fc2 . fr'N_i  can  be  ol  ained  itera¬ 

tively  as  follows. 


1.  Initial  values  of  k0  to  fr/v_j  re  chosen. 


2.  These  are  used  in  the  right  side  of  (15),  yielding 


Vo 

ko 

ki 

A'2 

k\- 2  ks~  i 

-1 

^[0] 

k[ 

1 

k  i 

ki 

h  . 

.  Ayv-i  0 

0V[1] 

k'7 

=  5(0) 

k2 

k'3 

*  4  . 

0  0 

0,[2] 

*V- 1  . 

A’.v-i 

0 

0  . 

t  d>?(.V  -  1]  _ 

3.  The  squared  distance  between  {A:,}  and  {A’(}  is  calculated. 

A’-i 

n=0 

4.  If  c  is  less  than  a  pre-established  tolerance,  k0,k\ . A\v_i  is  the  answer.  Otherwise,  let 

A'o  —  A*q,A.'i  =  k\, . . . ,  fc/y_]  =  k't- _ j ,  and  return  to  step  (2).  This  algorithm  is  not  guaranteed 
to  converge. 

Note  that  there  are  many  equations  that  can  be  employed  instead  of  (15). 

3.4.3  Design  of  the  power  spectrum 

Jn  this  section,  the  performance  of  the  SDMA  scheme  against  jamming  or  exploitation  attempts 
is  discussed.  Since  the  spectrum  of  the  SDMA  signal  is  dependent  on  the  wavelet  selected,  one 
question  arises  naturally:  “Docs  there  exist  a  wavelet  that  causes  the  spectrum  of  the  SDMA  signal 
to  look  like  that  of  white  noise?".  This  is  equivalent  to  the  requirement  that  <J>(w)  =  £2,n 
is  constant  within  a  fairly  wide  frequency  band.  This  requirement  is  vague  since  the  scalars  which 
m  selects  are  unknown  and  the  definition  of  “fairly  wide”  is  vague.  It  is  useful  to  discuss  a  related 
equation  that  gives  some  insights, 

S{ 2mw)  =  constant  (lb) 

Note  that  m  takes  on  all  integer  values  in  this  equation.  Defining  q(r)  =  S(2r),  (1G)  ran  be  further 
written  as 

J2  s^mu)  =  n  + m ) 

m€£  mtZ 


1G 


Lot 


G(z)  =  XI  y(x  +  m> 

m€?. 

Hence  equation  (16)  implies  G  is  constant  for  all  real  x.  Because  of  the  periodicity  of  G ,  it  can  be 


expanded  in  a  Fourier  scries, 


The  coefficients  are  given  by 


G(x)  =  £  cue2*" 
n€Z 


ri/i  ;i/2  ^ 

n=  G{x)t  2*>ni  dx=  ^S(i  +  m)e 

J- 1/2  J- 1/2  rr, 

-  L 


-  2jti ni 


»/*  mez 


— 2*»nx 


=  $(2*n.) 


where  3  is  the  Fourier  transform  of  <7.  Hence, 


G(*)  =  XIj(2*n)C2™* 

.,€/ 

The  condition  g(2rtm)  =  0  or  /K  S(2“')e~‘“2rm  dur  =  0  for  m  ^  0  is  necessary  and  sufficient  for  G  to 
be  constant,  hence  i6  necessary  for  S(‘2'"u>)  to  be  equal  to  a  constant. 


3.5  Bandlimited  orthogonal  wavelet  bases 


The  problem  of  channel  crosstalk  (intersymbol  interference)  that  arises  in  the  presence  of  asyn¬ 
chronous  transmitter  operation  and  propagation  delay  was  mentioned  above.  One  way  to  combat 
this  problem  is  to  place  different  transmitters’  signals  in  disjoint  frequency  bands,  ibis  section 
derives  -mill i- hand”  wavelet  bases  whose  elements  at  different  scales  occupy  disjoint  bands  in  the 


frequency  domain. 

One  such  basis  is  generated  by  the  (frequency  domain)  mother  wavelet 

(l  7T  <  |w|  <  2?r 

W(u)=  <  (17) 

I  0  otherwise 

With  a0  =  2  and  b0  =  1,  this  wavelet  is  easily  shown  to  generate  an  orthonormal  basis  of  /,2  [5].  This 
basis  gives  insight  about  how  to  construct  other  bandlimited  wavelet  bases  in  which  the  dilations  of 
the  mother  wavelet  W  do  not  overlap  in  the  frequency  domain. 

In  what  follows,  the  connected  components  of  the  support  of  W  will  be  called  “bands”.  J  he 
wavelet  defined  in  (1 7)  is  thus  a  “two- band”  wavelet.  The  following  sections  describe  the  construction 
of  “ri  band”  orthonormal  wavelet  bases. 
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3.5.1  Design  of  n-band  orthogonal  sets 

Suppose  the  support  So  of  \Y  is  the  union  of  n  disjoint  intervals 

So  =  f^Oi  ■S]]  ^3]  I^J  •  •  ■  l 1  ■Sn] 

Then  the  support  of  is 

Sm  =  [502-m^,2-'"j(J[522— ^32-m]|J...U[5T1_I2-m)5ri2-mj 

If  the  measure  of  Sm  flSm'  *s  zero  for  all  lion-equal  integers  m  and  m1,  Sq  is  called  an  orthogonal 
support,  lienee,  if  the  support  of  W  is  orthogonal,  then 

(W£,  l V»',)  =  ^(lV»,  <,)  =  0  for  all  m  /  m' . 

In  order  that  the  time-shifted  and  dilated  replicates  IV"  of  IV  form  an  orthonorma.l  set  in  L2, 
it  remains  to  assure  that  W™  and  IV"',  are  orthogonal  for  n/n'  and  all  m.  This  is  simplified  by 
the  observation  that 

{wz.w£)  =  {wz,ws') 

Thus  one  must  only  verify  that  the  time-shifted  replicates  of  IV  are  orthogonal  at  a  single  scale.  J  he 
following  two  theorems  describe  how  the  support  Sq  of  11 0  ran  be  selected  to  ensure  orthogonality. 

Theorem  3.1  Consider  real  numbers  0  <  Co  <  ci  <  c2  <  ■  •  •  <  c„_i  <  c„  =  2co  (or  0  >  Co  >  ej  > 
c2  >  ■■■  >  c„_,  >  c n  =  2 c0).  Then 

[c0i  2f0j  =  [co,ci](J[ci,C2](JSc2,C3](J  -  IJ[c„_i,2co] 

Let  pi , . . .  ,p„  e  and  dilate  each  sub-interval  [cj,cJ+i]  by  2p>+‘  to  form  a  set 

So  =  [c<j2p‘ ,  Ci2pl]  (J[ci2p?,c22p'](J  . .  - 1 2P” ,  2c02r”] 

=  (50,  «l]  (J^2’  53]  (J  •  -  -UN"’  S2’>-1) 

Then  So  is  orthogonal. 

Theorem  3.2  Suppose  a  support  So  is  orthogonal.  Then  Sq  can  be  generated  by  the  method  mentioned 
in  Theorem  3.1  when  parameters  are  properly  chosen  and  possibly  some  subintervals  are  deleted. 
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The  proof  of  these  theorems  are  given  in  (28j. 

Suppose  the  support  of  the  mother  wavelet  IV  is  orthogonal.  Then  {IV™}  will  become  an 
orthogonal  wavelet  set  if  and  only  if  IV  also  meets 


<W0\W0"')  =  0  for  n  ?  n  \  11,  n'  e  2 


(18) 


dt 


This  is  equivalent  to 

0  -  (Wo « wo) 

=  [  \V{t  -  n)W’(t  -  n') 

J  R 

=  ~  (  c-‘uw\V(u)c'n'uW'{uj)  du 

2 ?r  ./ r 

=  ~  (  e'("'-n^\\V(uj)\2 

2tr  Jk 

for  all  n  ^  n' .  With  k  —  n  -  n'  ^  0,  this  expression  becomes 

J  e^IVVfuOj2  du  =  0 

for  any  non-zero  integer  k.  Some  solutions  of  this  equation  are  discussed  below. 


(19) 


(20) 


3.5.2  Some  bases 

Suppose  VV  is  the  Fourier  transform  of  the  mother  wavelet  IV  and  its  support  on  the  frequency 
domain  is  symmetric  with  positive  part  is 


V  =  [c02Pl  ,  ci2PI]^J[c12P2,c22Pij(J  •  •  ■[J[(:r,-i‘£Pn *cn2Pn] 


where  0  <  cq  <  c\  <  c2  <  .  ■ .  <  c„  ~  2co  and  pi,p2,...,p„  are  integers.  Also  suppose  that,  the 
amplitude  of  VV  equals  1  in  V  and  that  the  phase  is  0.  Then  the  left  side  of  (20)  can  be  calculated 
as  follows: 


f  c'*"|VV(w)|2  du) 


=  1 

r-c02Pl 
{  1  + 

,Ct  }P1 
/  +  J 

j  + 

/  +■■ 

r-Cn_  jjsn  -Cl 

■  ■+  /  +  / 

VcoJSi  J 

/c,  2SS 

Jcn, 

rc\2^  rc2V’4  pen!''" 

=  2(/  +  /  +...+  /  )  cos  (ku)  du 


2 

=  jm{- 6in(c02Plfc)  +  sin(ci2Plfc)  -  sin(ci2PJfc)  ■+  sin(c22P3Jt)  - 

...  -  sin(ctl_,2p"*)  +  sin(crl2p"fc)}  =  0  (21) 
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i  C_  w 


This  equation  is  equivalent,  to  the  satisfaction  of  the  following  equation  for  any  non-zero  integer  k: 

-  sin(c02pi  k)  +  s'm(cx2p' k)  -  sin(ci2P2fc)  +  sin(c22P2/t)  - 

...  -  sin(cn_12r'1fc)  -f  sin(cn2p"fc)  =  0  (22) 

If  Co2pl,  ci2Pl,  C)2P2,  c22P2,  . . .,  ctl_i2Pn,  and  cn2Pn  satisfy  (22),  then  A/e02Pl,  Mcx 2P\  M cx 2P2, 
Afc22P2,  . . .,  Afcn_ \2r" ,  and  A/c„2Pn  will  also  satisfy  (22)  for  any  natural  number  A/.  So  there 
exists  a  support  of  a  “lowest.”  level,  which  will  be  called  the  mother  support.  In  other  words,  if  the 
end  points  of  the  support 

V  =  (c02p',c12p‘](J[ci2?,3,c22P2HJ...lJ[cn_12p",cn2p''] 
satisfy  (22)  and  the  end  points  of  the  support 

VM  =  (c02pl  / M,  c\2Pl  /M]  (J[c j2P2 /Af ,  c22Pj /jW  j  (J  . .  .[J[en_,2p"/^,  cn2*"/M) 
do  not  meet  (22)  for  any  natural  number  Af ,  then  V'  is  the  mother  support. 

Example  3.1  : 


Equation  (22)  is  difficult  to  solve  explicitly  because  it  is  nonlinear  in  several  variables.  It  car. 
be  simplified  by  considering  the  special  situation  where  the  support  consists  of  only  two  symmetric 
intervals;  i.e.,  the  support  is  [-2c02Pl ,  ~r02pi]  (J!co2Pl ,  2c02Pl]  (c0  >  0).  Denoting  a  =  c02Pl ,  (22) 
becomes 

-  sin(afc)  +  sin(2 ak)  =  0 


or 


sin(a.k)[2  cos(ak)  —  1]  =  0 


(23) 


The  only  valid  solution  of  (23)  is  a  =  mi  for  m  e  N  (the  other  solution  a  =  2mn  ±  ir/3  will  not 
satisfy  (23)  for  any  non-zero  integer  k ), 

Hence,  the  support  is  [-2mir,-mi]  U[mir,  2mir],  the  mother  support  is  [-2x,  -ir)  U(tt ,  2ir).  An 
orthogonal  wavelet  basis  can  be  constructed  whose  mother  wavelet  is 


!1  if  ir  <  |w|  <  2ir 
0  otherwise 


(2-1) 
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4  Research  Results:  Frequency  Hopping  Signal  Demodulation 

Research  in  this  area  was  focused  on  dernodulai ion  of  intercepted  frequency  hopping  spread  spec¬ 
trum  signals  in  an  uncluttered  or  lightly  cluttered  environment.  Several  slow  frequency  hopping 
signals  (digital  voice,  text,  analog  audio,  etc.)  were  generated  to  form  a  modest  library  of  test 
data.  This  collection  includes  both  analog  and  digital  data  and  various  suppressed  carrier  baseband 
modulation  techniques  are  represented.  Software  for  time  windowed  fast  Fourier  transform  (FFT) 
and  wavelet  transform  (WT)  processing  of  spread  spectrum  data  was  written.  A  simple  Viterbi- 
type  algorithm  that  operates  on  raw  time-frequency  (periodogram)  or  time-scale  (scaleogram)  data 
and  produces  tracks  of  frequency  hopping  signals  through  time-frequency  or  time-scale  space.  This 
software  was  used  to  support  empirical  evaluation  of  WT  measurements  as  the  basis  for  estimating 
the  spreading  sequence  of  a  single  slow  frequency  hopping  signal  in  a  noise-free  and  clutter-free 
environment.  Several  wavelets  were  tested,  and  it  became  evident  that  the  wavelets  that  provided 
scaleogram  measurements  that  led  to  the  best  spreading  sequence  estimates  were  essentially  sinu¬ 
soidal.  Indeed,  the  fact  that  frequency  hopping  signals  are  carried  by  sinusoidal  waveforms  seems  to 
make  Fourier  analysis  particularly  well  suited  to  this  problem.  However,  in  experiments,  a  hybrid 
algorithm  using  fixed- resolution  FFT  processing  and  multiresolution  WT  processing  provided  better 
spreading  sequence  estimates  than  cither  FFT  or  WT  processing  alone.  The  WT  component  of  the 
algorithm  provided  better  localization  of  the  hopping  times  of  the  signal  than  the  fixed  resolution 
FFT  component  (writh  relatively  long  integration  period  chosen  to  provide  accurate  frequency  mea¬ 
surements)  than  either  WT  or  FFT  processing  alone.  In  addition,  the  efficiency  of  fast  algorithms 
for  computation  of  wavelet  transform  values  at  dyadically  arranged  points  in  the  phase  plane  pro¬ 
vided  a  fast  approach  to  "focusing  in”  on  the  true  hopping  times  and  frequencies.  Some  additional 
experiments  with  lightly  cluttered  communication  environments  were  undertaken. 

Among  the  main  reasons  for  the  use  of  frequency  hopping  spread  spectrum  communication 
systems  is  their  robustness  with  respect  to  detection  and  jamming  [20],  Typical  slow  frequency 
hopping  signals  consist  of  an  information  signal  of  relatively  small  bandwidth  modulating  a  carrier 
signal  having  piccewisc-constant  frequency  that  varies  (“hops”)  within  a  band  much  larger  than 
that  of  the  information  signal  (figure  5).  The  spreading  function  is  known  to  intended  receivers. 


22 


Frequency  (KHz) 


1400 

noo  - 


I _ L. _ 1 _ i  .  .  r 

100  200 

Time  (ms) 

Figure  5:  A  typical  spreading  function  for  a  slow  frequency  hopping  communication  signal.  In  this 
context,  the  adjective  “slow”  means  the  period  between  hops  (which  may  be  fixed  or  variable)  is 
much  longer  than  the  duration  of  a  symbol  representing  a  single  binary  digit. 

but  unknown  to  unintended  receivers  who  may  wish  to  demodulate  or  jam  the  transmission.  This 
paper  presents  preliminary  results  of  research  into  techniques  for  estimating  the  spreading  function 
of  an  unknown  frequency  hopping  transmission  using  only  received  (intercepted)  data.  Since  such 
functions  arc  characterized  by  the  times  and  frequencies  associated  with  the  hops,  techniques  are 
sought  to  accurately  estimate  these  parameters.  Good  estimates  of  the  spreading  function  are 
directly  useful  in  demodulation  of  such  signals  and  may  also  be  useful  in  jamming  them. 

4.1  Time-frequency  techniques 

The  analysis  tools  used  in  this  work  were  the  short-time  Fourier  transform  (STFT)  and  the  wavelet 
transform  (WT).  Other  approaches  for  analysis  of  signals  with  time- varying  characteristics,  such  as 
the  Wigner-Villc  and  related  time-frequency  distributions  [9]  may  also  be  applicable  to  this  problem. 
These  were  not  used  in  the  work  presented  here,  however. 

4.1.1  The  STFT 

As  a  signal  analysis  tool,  the  primary  utility  of  the  Fourier  transform  is  decomposition  of  a  finite- 
energy  signal  into  individual  frequency  components.  Its  value  in  analysis  of  signals  with  time 
varying  characteristics,  such  as  frequency  hopping  signals,  is  limited  because  information  about 
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when  the  various  frequency  components  are  present  in  a  signal  is  embedded  in  the  phase  of  its 
Fourier  transform  in  a  way  that  is  not  easy  to  interpret.  The  magnitude  of  the  Fourier  transform 
of  a  frequency  hopping  signal,  for  example,  contains  no  information  about  when  the  hops  occur. 

The  STFT  or  time-windowed  Fourier  transform  of  a  finite-energy  signal  /:"<—♦  "A  is  a  function 
F  ;  K2  — *  E  defined  by 

F{t,u>)  =  J  /(r)te(r  -  t)e~'WT  dr 

In  this  expression,  w  is  a  bounded  function,  generally  even  and  with  bounded  support,  called 
a  window.  In  general,  window  selection  is  based  on  numerous  criteria  [15].  In  the  context  of 
frequency  hopping  signal  analysis,  the  primary  tradeoff  is  in  selecting  the  time  duration  of  the 
window.  Windows  of  short  duration  allow  precise  temporal  localization  of  hops  at  the  expense 
of  frequency  resolution.  Analysis  using  long  windows  sacrifices  temporal  resolution  but  improves 
frequency  resolution. 

4.1.2  The  WT 

The  wavelet,  transform  of  a  continuous-time  finite  energy  signal  /  with  respect  to  an  analyzing 
wavelet  y  has  values 

In  this  expression,  g‘  denotes  the  complex  conjugate  of  g,  n  >  0  is  called  the  dilation  or  scale 
parameter,  and  b  6  'S.  is  the  time  shift  parameter. 

The  wavelet  transform  has  been  shown  to  be  of  particular  value  in  detecting  and  temporally 
localizing  discontinuities  and  other  sudden  changes  in  signals  [1,  22].  The  ability  of  the  WT  using 
any  given  wavelet  to  resolve  narrowband  signals  that  are  closely  spaced  in  frequency  is  better  at  low- 
frequencies  than  higher  frequencies.  This  characteristic  is  ideally  suited  for  certain  applications  in 
which  frequency  resolution  requirements  diminish  in  proportion  to  the  center  frequency  of  the  band 
being  analyzed  (e.g.,  “constant-Q”  settings).  It  is  not  desirable  for  analysis  of  frequency  hopping 
signals,  however,  because  the  need  for  frequency  resolution  is  nnilorm  across  a  wide  bandwidth. 
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Figure  6:  Flow  diagram  of  a  hybrid  algorithm  for  estimating  t  he  spreading  function  of  an  intercepted 
frequency  hopping  bignal. 


4.2  An  estimation  algorithm 

The  characteristics  just  discussed  suggest  a  hybrid  scheme  for  spreading  function  estimation  in 
which  WT  analysis  is  used  to  estimate  the  hopping  times  and  STFT  analysis  is  used  to  estimate 
the  frequencies  of  the  signal  segments.  Such  an  algorithm  is  depicted  in  figure  6. 

At  the  first  stage,  the  WT  of  a  segment  of  the  intercepted  data  is  computed  digitally.  The 
discrete  wavelet  transform  used  at  this  stage  is  dyadic  and  provides  one  output  value  mo.j  at  the 
coarsest  scale  a o-  It  provides  two  output  values  w j,i  and  w]  2  at  the  next  finer  scale  aj  =  2_,ao 
At  scale  an  =  2_ntioi  it  yields  2n  values  tnnj,  resulting  in  a  dyadic  tree  of  output  values 

ending  at  the  finest  scale  a ^  =  2~noq.  The  occurrence  of  a  hop  wit  hin  the  segment  is  determined  by 
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comparing  the  sequences  |u'o.i),  and  | w’o.i | ,  |ti-i  2|,  ...,  |ie/V 2*|.  These  sequences  may 

be  interpreted  as  crude  local  spectral  estimates  of  the  signal  near  the  leading  and  trailing  edges 
of  the  segment,  respectively.  Thus,  if  a  hop  occurs  within  the  segment,  the  change  in  frequency 
structure  of  the  signal  will  be  reflected  in  the  sequences.  If  a  hop  is  indicated,  the  outputs  form  the 
finer  scales  are  used  to  estimate  the  precise  lime(s)  of  the  hop(s). 

The  next  stage  of  the  algorithm  processes  a  segment  of  data  bounded  in  time  by  two  hops,  as 
detected  in  the  first  stage  of  the  algorithm,  using  a  ST1T  with  a  window  whose  length  precisely 
matches  that  of  the  segment.  This  is  implemented  by  a  zero-padded  FFT  [18].  The  spectral  estimate 
obtained  is  used  to  estimate  the  frequency  of  the  carrier  between  hops. 

4.3  Simulation  results 

A  digitized  frequency  hopping  signal  with  the  following  characteristics  was  produced: 

•  BPSK  baseband  modulation 

•  Bit  rate  of  9600  bps 

•  Hopping  rate  varying  randomly  between  50  and  250  ms 

•  Smallest  frequency  hop  of  20  Kllz 

•  No  noise  or  interfering  signals 

Because  the  true  spreading  function  of  this  signal  was  known,  it  was  possible  to  reconstruct  the 
baseband  signal  from  the  frequency  hopping  signal  simulating  the  reconstruction  possible  under  ideal 
conditions  at  an  intended  receiver.  The  signal-to-noise  ratio  (SNR)  obtained  in  this  reconstruction 
was  approximately  25  dB. 

The  hybrid  algorithm  described  above  was  used  to  estimate  the  spreading  function  of  the  sim 
ulated  signal.  The  parameters  of  the  algorithm  were  chosen  based  on  the  assumption  that  the 
shortest  time  between  hops  would  be  no  less  than  25ms  and  that  the  smallest  frequency  hop  would 
be  no  less  than  10  KHz.  Otherwise,  the  algorithm  used  no  specific  a  prion  knowledge  about  the 
signal.  The  baseband  signal  was  reconstructed  using  the  estimated  spreading  function  with  a  SNR 

of  approximately  14  dB. 
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A  similar  spreading  function  estimation  algorithm  using  STFT  rather  than  WT  to  estimate 
the  hopping  times  was  implemented  and  run  on  the  same  data.  It  was  found  to  be  much  more 
computationally  eflicient,  but  the  csti  .ated  hopping  times  were  not  as  accurate  as  those  obtained 
using  WT.  T  his  loss  of  accuracy  reduced  the  reconstruction  SNR  to  approximately  8  dB. 

4.4  Demodulation  versus  jamming 

In  principle,  accurate  estimation  of  the  spreading  function  of  a  frequency  hopping  signal  can  be 
useful  in  jamming  as  well  as  in  demodulation.  For  jamming  applications,  a  Imp  must  be  detected 
and  the  frequency  of  the  new  segment  estimated  essentially  immediately  in  order  for  the  jamming 
transmitter  to  switch  to  the  new  frequency  before  the  next  hop.  Thus,  estimation  algorithms  must 
be  causal  and  very  fast  for  this  application.  On  the  other  Hand,  the  accuracy  of  the  estimate  may 
not  be  crucial.  Just  having  a  general  idea  where  in  the  phase  plane  the  signal  is  concentrated  will 
allow  much  more  efficient  use  of  jammer  energy  than  constant  jamming  of  a  broad  frequency  band. 

For  demodulation  of  an  intercepted  signal  neither  causality  nor  fast  algorithm  execution  may  be 
necessary.  Latencies  of  seconds  or  even  minutes  may  have  little  effect  on  the  utility  of  the  spreading 
function  estimate  in  demodulation.  However,  accuracy  of  the  estimate  in  this  application  is  crucial. 

5  Publications 

Kesults  from  this  research  effort  have  been  published  as  follows. 

5.1  Journal  and  conference  proceedings  papers 

J  he  following  paper'  describing  results  obtained  under  this  effort  were  published  or  submitted  for 
review  and  possible  publication: 

1.  N.  Bhouti  and  D.  Cochran,  “Multiscale  Time-Frequency  Techniques  for  Spread  Spectrum  De¬ 
modulation  and  Jamming,”  Proceedings  of  the  lu-enty-  Sixth  Asilomar  Conference  on  Signals, 
Systems,  and  Computers ,  October  1992. 

2.  1).  Cochran  and  C.  Wei.  “Scale  Based  Coding  of  Communication  Signals,”  Proceedings  of 
the  IFFF-SP  International  Symposium  on  Time- Frequency  and  Time-Scale  Analysis,  October 
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X  ])  Cochran  and  S.  Enscrink,  “Use  of  Cvclostni  ion  ary  Features  in  Frequency  Hopping  Signal 
Segment  Association,”  Proceedings  of  t he  Workshop  on  Cyclostationary  Signals,  August  1992. 

4.  C.  Wei  and  I).  Cochran,  “Construction  of  Discrete  Orthogonal  Wavelet  Bases,”  Record  of  the 
IEEE  International  Symposium  on  Information  Theory,  January  1993. 

5.  !).  Cochran  and  C.  Wei,  “Scale-Division  Multiple  Access  -  Part  I:  Wavelet  Coding  and  Spread 
Spectrum,”  Submitted  to  IEEE  Transactions  on  Information  Theory. 

5.2  Theses 

The  following  thesis  describing  results  obtained  under  this  effort  was  prepared  and  defended: 

•  N.  Uhouri,  Application  of  Multircsolution  Time- Err. qucncy  Techniques  to  Spread  Spectrum  De¬ 
modulation  and  Jamming,  M.S.  Thesis,  Arizona  Slate  University,  December  1991. 

An  additional  thesis  is  in  preparation  and  is  expected  to  be  defended  in  May,  1993: 

•  C.  Wei,  Scale-Division  Multiple  Access,  M.S.  Thesis,  Arizona  State  University,  May  1993  (in 
preparation). 

6  Participating  Personnel 

Personnel  contributing  to  this  research  effort  were: 

1.  D.  Cochran,  Principal  Investigator 

2.  C.  Wei,  Research  Assistant 

3.  N.  Uhouri,  Research  Assistant 

4.  D.  Sinno,  Research  Assistant 

The  research  component  of  the  following  graduate  degree  was  undertaken  as  part  of  this  effort: 

l.  N.  Bhouri,  Master  of  Science  in  Electrical  Engineering,  December  1991.  Thesis:  Application  of 
Muitiresoluiion  Time-Erequency  Techniques  to  Spread  Spectrum  Demodulation  and  Jamming 


2-  Wei,  Master  of  Science  in  Electrical  F.ngi nceri  ng , ToTom plot o  May  1093.  Thesis:  Sralc- 
Division  Multipit  Access 

7  Research  Interactions  • 

The  following  research  interactions  were  related  to  this  project: 

1.  Alt  oral  presentation  entitled  “Application  of  Multiscale  Processing  to  Characterize  Frequency 
Hopping  Spread  Spectrum  Signals”  was  prepared  by  D.  Cochran  and  N.  Bhouri  presented  by 
D.  Cochran  at  the  AFOSR/AFIT  Symposium  on  Application  of  Wavelets  to  Signal  Processing 
at  Wright-Patterson  AFB  in  March  1991. 

2.  D.  Cochran  attended  a  series  of  meetings  at  Wright  Laboratories  in  October  1991  to  dis 
cuss  this  project  and  related  topics  with  the  AFOSR  Program  Manager  (J.  Sjogren),  An 
Force  Laboratory  personnel  (.1.  Stephens,  L.  Gutman,  J.  Tsui),  and  AF11  Faculty  (M.  Oxley 
G.  Warhola).  In  particular,  a  progress  briefing  on  this  project  was  presented  to  the  Program 

Manager  at  this  time. 

"  D.  Cochran  presented  an  invited  lecture  in  the  AFJT/ AFOSR  Distinguished  Lecture  Serin 
at  the  Air  Force  Institute  of  Tpthnologv  in  September  1992.  During  this  visit  to  Wright 
Patterson  AFB,  he  collaborated  extensively  with  AFIT  faculty  (M.  Oxley,  !L  Suter)  and  me 
with  Wright  Laboratories  personnel. 

•1.  1).  Cochran  presented  talks  on  this  project  at  (a)  the  University  of  Arizona  lectrical  Engined  - 
ing  Department  Colloquium  [35  attendees)  in  Novemli£l_lR90,  (b)  the  IEF  .  Signal  Processing 
and  Communications  Society  Phoenix  Chapter  Meeting  [48  attendees)  in  February  1991,  (c. 
the  Motorola  Government  Electronics  Group  Signal  Processing  Seminar  Series  [13  attendees) 
in  April  1991,  and  (d)  MIT  Research  Laboratory  for  Electronics  Signal  Processing  Group 
Seminar  [12  attendees]  in  February  1992. 

8  Inventions  and  Patents 

No  inventions  or  patent  disclosures  resulted  from  this  research  program. 
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9  Additional  Information 

Work  on  the  most  promising  aspects  of  this  project  is  continuing  under  grant  No.  P4  19620-93-1-005 
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